function inner=obj(coeff,P1,P2,epsilon_sim)
    comp1 = [P1 P2]*coeff;
    comp2 = mean(max(epsilon_sim+[coeff.' 0],[],2));
    inner = -(comp1-comp2);
end